Thermal conductivity of pure and Zn-doped LiCu 2 02 single crystals 
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We report a study of the low-temperature thermal conductivity (k) of pure and Zn-doped LiCu2 02 
single crystals. The k(T) of pure LiCu2C>2 single crystal shows a double-peak behavior, with two 
peaks locating at 48 K and 14 K, respectively. The different dependences of the peaks on the Zn 
concentration indicate that the high-T peak is likely due to the phonon transport while the low-T 
one is attributed to the magnon transport in the spin spiral ordering state. In addition, the magnetic 
field can gradually suppress the low-T peak but does not affect the high-T one; this further confirms 
that the low-T peak is originated from the magnon heat transport. 

PACS numbers: 66.70.-f, 75.47.-m, 75.50.-y 
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I. INTRODUCTION 

LiCu202 is the first example of Cu-based multiferroic 
material and is particularly attractive because of its one- 
dimensional spin structured— It promises a new routine 
to find multiferroicity in some low-dimensional quantum 
magnets that exhibit the magnetic frustration effect 
It is known that LiCu202 crystallizes in an orthorhombic 
unit cell with space group Pnma and lattice parameters a 
= 5.734(4) A, b = 2.856(2) A and c = 12.415(6) A at room 
temperature r 1 There are an equal number of Cu + and 
Cu 2+ ions in distinctly nonequivalent crystallographic 
positions, only the latter of which carry spin S = 1/2. 
The Cu 2+ ions are sitting on the center of edge-sharing 
Cu04 plaquettes and form edge-shared chains running 
along the b axis with the Cu-O-Cu bond angle of 94° #1 
The competition of the nearest-neighboring ferromag- 
netic (FM) interaction and the next-nearest-neighboring 
antifcrromagnetic (AF) interaction of Cu 2+ spins leads 
to magnetic frustration and a spiral (helicoidal) magnetic 
order below ~ 24 K.— More exactly, two AF transitions 
were found at T NX = 24.6 K and T N2 = 23.2 K with 
a sinusoidal spin order at T/y 2 < T < Tjyi and an in- 
commensurate cycloidal spin order at T < T/V2<^r— A 
spontaneous polarization along the c axis emerges at the 
second phase transition^ and was discussed to be orig- 
inated from the inverse DM interaction of neighboring 
spins or the nonrelativistic exchange. 3,7 Furthermore, the 
magnetic field applied along the b axis leads to the Cu 2+ 
spin spiral plane flipping from the be to the ab plane and 
consequently results in the flip of the polarization from 
the c to the a direction. 1 

The low-temperature thermal conductivity is an ef- 
fective probe for studying the transport properties of 
phonons and magnetic excitations, which might be of 
particularly interesting for the low-dimensional quan- 
tum magnets^ — It is well known that the elemen- 
tary excitations in magnetic long-range ordered states, 
magnons, can contribute to the heat transport properties 
by acting as either heat carriers or phonon scatterersJ^ 



In low-dimensional spin systems, even without exhibit- 
ing long-range ordering, the magnetic excitations can ef- 
fectively transport heat because of the strong quantum 
fluctuations*^ For example, the recent experiments re- 
vealed extremely large magnetic heat conductivity in one- 
dimensional spin 1/2 systems, such as SrCuC>2, S^CuOs, 
CaCu 2 3 and (Sr,Ca,La) 14(^240-41^— These results 
were theoretically well understood as the ballistic trans- 
port of spinons or magnons. On the other hand, the 
thermal conductivity can effectively detect the transi- 
tions of magnetic structure like spin flop, reorientation or 
polarization^ In particular, the low-T thermal conduc- 
tivity was found to display drastic changes across these 
kinds of transition in some other families of multiferroic 
materials, for example, HoMn03 and GdFe03i 25 ' 26 

Although the low dimensionality of the magnetic struc- 
ture LiCu202 has already known for long time, the heat 
transport properties have not been investigated for this 
material. In this work, we study the temperature and 
magnetic-field dependences of thermal conductivity (k) 
of LiCu202 single crystals for probing the nature of mag- 
netic excitations and the coupling between spin and lat- 
tice. It is found that the k(T) data show two peaks at 48 
K and 14 K, which are above and below the long-range 
magnetic transition temperatures, respectively. The Zn 
substitution for Cu is found to be able to effectively sup- 
press the two peaks but show quite different doping de- 
pendences. The magnetic field applied in the ab plane 
only suppresses the low-T peak. These results manifest 
that the low-T peak is likely due to the magnons con- 
tribution to the heat transport acting as heat carriers, 
while the high-T peak is phonon peak. 



II. EXPERIMENTS 

High-quality single crystals of LiCu2- 2: Zn 2 ;02 with the 
nominal compositions x = 0, 0.02, 0.04, 0.10 and 0.20 are 
grown by a self-flux method. Correspondingly, the actual 
Zn contents of these crystals are x = 0, 0.013, 0.027, 0.064 
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FIG. 1: (Color online) (a) X-ray rocking curve of the (006) 
reflection of a LiCu202 single crystal. The FWHM of the 
peak is only 0.05°. (b) X-ray (00Z) diffraction pattern of the 
same crystal. Inset: The photograph of this crystal. 

and 0.177, measured by the inductively-coupled plasma 
atomic-emission spectroscopy (ICP-AES) (with 10% un- 
certainty of measurements). The as-grown LiCu2 02 sin- 
gle crystals have plate-like shape with typical size as large 
as 10 x 8 x 0.5 mm 3 . Upon doping Zn, the sizes of crys- 
tals decrease gradually to about 5 x 3 x 0.5 mm 3 for x 
= 0.177. The quality of the crystals, judged from the x- 
ray diffraction results, does not decay significantly. The 
largest surfaces of single crystals are confirmed to be the 
ab planes by the x-ray diffraction and Laue back reflec- 
tion. 

The magnetic susceptibility (x) is measured using a su- 
perconducting quantum interference device magnetome- 
ter (SQUID, Quantum Design). The specific heat is mea- 
sured by the relaxation method in a physical property 
measurement system (PPMS, Quantum Design) from 2 
to 300 K. The in-plane thermal conductivity is measured 
in PPMS with "one heater, two thermometers" configu- 
ration or in a 4 He cryostat using a Chromel-Constantan 
thermocouple^ - — 

III. RESULTS AND DISCUSSION 

LiCu2 02 single crystals have plate-like shape and shin- 
ning surfaces. Using x-ray diffraction, it is found that the 
large surface of the crystals are the ab plane, so it is easy 
to get the (00?) diffraction pattern. Figure 1 shows the x- 
ray diffraction results of a representative LiCu2 02 single 
crystal, including the (00/) diffraction pattern and the 
rocking curve of (006) peak. Figure 1(b) does not show 
any sign of other phases, indicating the high purity of this 
sample. The full width at half maximum (FWHM) of the 
(006) reflection is as small as 0.05°, indicating the perfect 
crystallinity of this sample. X-ray Laue back reflection is 
also used to determine the crystallographic axes. A fine 
twin structure of the ab plane is found, which is in agree- 
ment with the previous observation through a polarized 
optical microscope^ The origin is that the a-axis lattice 
constant is nearly the twice of the 6-axis length. Because 
of the existence of twin structure, it is impossible to mea- 
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FIG. 2: (Color online) (a,b) Temperature dependences of 
dx/dT for LiCu2- a: Zn a; 02 single crystals along the ab plane 
and the c axis measured with 1000 Oe magnetic field, (c) 
Low-temperature specific heat of LiCu2-iZn I 02 single crys- 
tals. For clarity, the data are shifted upward by 2 J/Kmol 
one by one. 



sure the in-plane anisotropy of the transport properties 
of LiCu202 single crystals, although all the samples are 
cut along the a (or b) axis. 

The magnetic susceptibility and specific heat are mea- 
sured to characterize our crystals. It has been known 
that the main feature of x(^) f° r pure LiCri202 is a 
broad maximum at ~ 36 K, which is a characteristic of 
quasi-one-dimensional magnet, and the onsets of long- 
range magnetic orders at T/vi and TV 2 induce a sharp 
anomaly at the temperature derivative of the magnetic 
susceptibility, dx(T)/dT, for magnetic field along the c 
axis and the ab plane, respectively— These results arc 
well reproduced in our crystals, as shown in Figs. 2(a) 
and 2(b). Furthermore, upon doping Zn, the peak posi- 
tions of dx{T) I dT gradually move to lower temperatures, 
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FIG. 3: (Color online) Thermal conductivity of pure and Zn- 
doped LiCu2C>2 single crystals. 

indicating the suppression of long-range magnetic order 
with increasing nonmagnetic impuritiesi 10 ! 11 Figure 2(c) 
shows the low-temperature specific heat of pure and Zn- 
doped LiCu2- a; Zn 2: 02 single crystals. For pure sample, 
the C(T) curve shows two small but clear peaks at Tvi = 
24.7 K and T/V2 = 23.2 K, respectively, which are known 
to be due to the magnetic phase transitions from para- 
magnetic state to a sinusoidal spin ordering and then to 
a helicoidal spin ordering.^— Upon doping Zn, these two 
peaks shift to lower temperatures and become weaker 
and their positions have good correspondence with those 
in dx{T)/dT. When x arrives 0.177, the two peaks in 
C(T) are not distinguishable from each other and evo- 
lute to a hump-like anomaly at ~ 17 K. This evolution 
of the specific heat with Zn doping is compatible with 
some earlier reports ) 10 ^ 1 in which it was discussed that 
Zn doping results in either the phase transition of short 
range or significant inhomogeneity. 

Figure 3 shows the low-temperature thermal conduc- 
tivity of pure and Zn-doped LiCu2 02 single crystals. The 
thermal conductivity of pure crystal is rather large but 
its temperature dependence is apparently different from 
that of usual insulators. 18 With lowering temperature, 
the k increases quickly and achieves a high value of 57 
W/Km at 48 K. A remarkable feature of n(T) is the ap- 
pearance of two peaks at 48 K and 14 K. It is notable 
that the minimum between two peaks is located at ~ 
22 K, having a good correspondence to the positions of 
specific-heat peaks. This suggests a possible origin of 
the double peaks of n(T), that is, their appearance is 
due to the formation of a minimum, caused by strong 
phonon scattering by the spin fluctuations at the criti- 
cal regions of magnetic transitions.— This phenomenon 
was also found in some multiferroic materials exhibit- 
ing strong spin-phonon coupling, such as HoMn03<2ti Al- 
though this explanation does not bring obvious contra- 



dict with the impurity (Zn) doping effect, however, it 
meets with difficulty in understanding the magnetic-field 
dependence of thermal conductivity, as shown below. 

It is found that the two k(T) peaks show rather dif- 
ferent dependences on the concentration of Zn. While 
the high-T peak weakens gradually with increasing x, the 
low-T one shows a weak x dependence for low doping lev- 
els x < 0.027 but it is so strongly damped for x > 0.064 
that it is almost smeared out completely. Note that the 
evolution of low-T peak seems to have some direct rela- 
tionship to the Zn-doping effect on the specific-hcat data. 
As can be seen in Fig. 2, the two specific- heat peaks are 
also smeared out for x > 0.064, which manifests that the 
long-range spin orderings are destroyed by such high im- 
purity dopings. This comparison naturally suggests that 
the low-T peak of t(T), which locates at temperature 
just below the phase transition of long-range spin order- 
ing, is likely due to the heat transport by the magnons 
in the antiferromagnetically ordered state. 

As far as the high-T peak of k(T) is concerned, one pos- 
sible origin is the phonon peak as in usual insulators. 18 
The magnitude of phonon peak is known to be strongly 
dependent on the impurities and point defects in crystals. 
Apparently, the strong impurity dependence of high-T 
peak of LiCu2C>2 is compatible with such standard be- 
havior. However, one may note that the position of this 
"phonon peak" is located at a bit too high tempera- 
ture, compared to many other materials which presenting 
phonon peak below 20 K*i& Therefore, we need to con- 
sider another possible origin of the high-T peak due to 
the magnetic excitations transporting heat, as evidenced 
in many low-dimensional spin systems, such as SrCuC>2, 
Sr 2 Cu0 3 , CaCu 2 3 , Sr 14 Cu 2 40 4 i, LasCuO^-S^^ 
In these materials, the magnetic term of thermal con- 
ductivity can be much larger than the phononic term 
and is also sensitive to the impurities. In particular, the 
magnetic excitations heat transport in quasi-one dimen- 
sional spin-1/2 materials was predicted to be the ballistic 
type. One direct experimental evidence for the ballistic 
transport was the non-magnetic impurity doping effect, 
in which the mean free path of magnetic excitations was 
found to be very close to the average distance between 
spin defects^ In this regard, the magnetic heat transport 
in these low-dimensional spin systems is usually obtained 
from the strong anisotropy of thermal conductivity j 16 i 17 
that is, one can get the purely magnetic thermal con- 
ductivity by subtracting the phonon term, which is es- 
timated from the thermal conductivity perpendicular to 
the one-dimensional spin chain or the two-dimensional 
spin network. It will be interesting to get the magnetic 
heat transport in LiCu202 by using the same method and 
compare the Zn-doping effect with that in other spin-1/2 
chain compounds. However, because of the in-plane twin 
structures of LiCu202 crystals, one cannot separate the 
heat transport along the a and the b axes. In addition, 
the heat transport along the c axis is also found to be im- 
possible to measure reliably because of the easy cleavage 
of the LiCu2C>2 crystals along the ab plane. It is therefore 
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FIG. 4: (Color online) (a) Temperature dependences of ther- 
mal conductivity of LiCu2 02 single crystal in 0-14 T. The 
direction of magnetic field is along that of the heat current, 
(b) Low-temperature k(H) isotherms at 5 K and 18 K. 



not feasible for us now to get the magnetic heat trans- 
port along the spin chains in these LiCu202 crystals. An 
effective way to detwin the LiCu2C>2 crystals is called 
for the investigation on the possible role of the magnetic 
excitations transporting heat. 

The effect of magnetic field on the thermal conductiv- 
ity is studied for in-plane fields up to 14 T. First of all, as 
shown in Fig. 4(a), the high-T peak is completely inde- 
pendent of the magnetic field, which would not be unrea- 
sonable if the high-T peak is a pure phononic behavior. 
For some well-studied low-dimensional magnetic materi- 
als, it is also found in recent experiments that the mag- 
netic heat transport is insensitive to the external mag- 
netic field, which is however due to the large exchange 
coupling in these materials, typically being of the order 
of magnitude of 100 meVi 22 ' 32 The exchange coupling in 
LiCu202 is known to be more than one order of mag- 
nitude smaller— and is therefore not much larger than 
the energy caused by the magnetic field in order of 10 
T. In this sense, the insensitivity of the high-T thermal 
conductivity to the magnetic field may not support the 
conjecture that the magnetic excitations are responsible 



for transporting heat above the magnetic phase transi- 
tions and the formation of high-T peak. Second, the 
"dip" between two k(T) peaks does not show any change 
in applied magnetic field, which immediately rule out 
the possibility that the "dip" is caused by strong spin- 
phonon scattering. In contrast, it has already been found 
in many AF materials that the strong magnetic field can 
suppress the spin fluctuations and recover the thermal 
conductivity if the spin-phonon coupling is considerably 
strong i 2 ^ 

On the other hand, the low-T peak is gradually sup- 
pressed by the magnetic field, similar to those in many 
antifcrromagnetically ordered materials^ The detailed 
field dependences of k are shown in Fig. 4(b), in which 
two k(H) isotherms at 5 and 18 K are included. It is clear 
that the thermal conductivity is monotonically decreased 
with increasing field, without showing any signature of 
saturation or drastic transition up to 14 T. This kind of 
field dependence is expectable for the magnon heat trans- 
port since the magnons tend to be less populated with 
increasing magnetic field. In addition, the field depen- 
dence at 18 K, which is near the low-T peak of k(T), is 
stronger than that at 5 K. It is also understandable be- 
cause the magnetic heat conductivity is apparently much 
larger at 18 K. This result shows a clear evidence that 
below the long-range-order transition temperature the 
magnons act as heat carriers in LiCu202- Note that both 
the k(T) and k(H) behaviors indicate that the magnon 
heat transport seems to be weakened with lowering tem- 
perature. This is mainly due to the decrease of magnon 
population and is reasonable for LiCu2 02 since there is a 
1.4 meV gap in the magnetic excitation spectrum, found 
by the electron spin resonance experiments^ It can be 
seen that the heat transport properties of LiCu202 are 
rather conventional^ without showing any peculiar be- 
havior of the low-dimensional quantum magnets. The 
reason is likely related to the rather strong spin frustra- 
tion in this material. 

It has been known that the in-plane magnetic field can 
rotate the spin directions and produce some spin-flop- 
like transitions at 2 Ti However, the k(H) curves do not 
show any anomaly across these transitions, in contrast to 
some observations in other compounds! 2 ^— One reason 
for the drastic change of k at the spin-flop transition is 
that the magnons are significantly populated because of 
the closure of the anisotropy gap^ Therefore, it is not 
clear whether the magnon spectrum is gapless at the spin 
re-orientation transition in this spirally ordered antifer- 
romagnet. 



IV. SUMMARY 

The Zn-doping and magnetic-field dependences of ther- 
mal conductivity of LiCu202 single crystals are care- 
fully studied. The k(T) data show a double-peak phe- 
nomenon. The higher-T peak at 48 K is due to the 
phonon heat transport, while the lower-T peak at 14 K 
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is a result of magnon heat transport showing up in the 
magnetic long-range-ordered state. The present results 
indicate that the magnetic heat transport of LiCu2 02 
behaves similarly as that in the three-dimensional anti- 
ferromagnets. The absence of the characteristic trans- 
port properties of that in the low-dimensional quantum 
spin systems may be related to the complexity of spin 
structure and spin frustration. 
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